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1.0  IDENTIFICATION 

1.1  Title 

Line  Printer  Plots  of  Two-Dimensional  Surfaces 

1 . 2  Identification  Name 
J5-NRL-SURFPL0T 

1 . 3  Classification  Code 
J5-0utput,  Plotting 

1 . 4  RCC  Identification  Number 
J  5008000 

1 • 5  Entry  Points 
SURFPL0T 

1 . 6  Programming  Language 

Language:  3600/3800  FORTRAN 

Routine  Type:  Subroutine 
Operating  System:  DRUM  SC</PE  2.1 

1 . 7  Computer  and  Configuration 
CDC-3800 

1 . 8  Contributor  or  Programmer 

Michael  A.  Tamny,  Code  5365T,  Airborne  Radar  Branch, 
Radar  Division 

1 . 9  Contributing  Organization 

NRL  -  Naval  Research  Laboratory, 

Washington,  D.C.  20375 

Manuscript  submitted  February  15,  1974. 
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1.10  Program  Availability 

1.10.1  Submittal : 

Program  write-up,  Fortran  source  deck 
with  comments,  source  listing. 

1.10.2  On  File: 

RCC  Program  Library 

1.11  Verification 

A  sequence  of  Cassini  Ovals,  defined  by  the 

equation  [  (x  +  a)  ^  +  y^j[(x  -  a)^  +  y^]  =  k\ 
where  k  equals  the  amplitude  and  a  equals  the 
constant  3.0,  was  used  to  test  SURFPL0T.  Several 
different  subroutine  parameter  values  were  used. 
The  resultant  plots  compared  favorably  with  those 
produced  by  the  contour  plotting  program 
J5-0RNL-C0NT0UR  on  the  CalComp  plotter. 

1.12  Date 
November  1973 

2.0  PURPOSE 

2 . 1  Description  of  the  Routine 

The  routine  SURFPL0T  generates  a  plot  of  two- 
dimensional  surfaces  using  the  line  printer. 

The  plot  is  a  plan  view  of  the  surface;  the 
amplitude  at  each  (x ,  y)  point  is  represented 
by  a  shade  of  gray.  The  gray  shading  is  created 
by  overprinting  with  the  line  printer  using  three 
line  printer  strikes  for  each  finished  line  in  the 
plot.  The  cost  in  turn-around  time  using  this 
routine  is  modest  as  it  can  use  either  of  the  NRL 
CDC-3800  systems  and  it  avoids  the  queue  at  the 
CalComp  plotter.  There  is  logic  to  handle  arrays 
that  are  too  large  to  or  do  not  reside  in  memory 
at  one  time  by  doing  the  plotting  piecemeal.  The 
plotting  is  ordinarily  accomplished  with  minimal 
user  involvement,  however  the  user  can  control 
the  nature  of  the  plot.  The  array  of  f(x,y) 
values  plotted  is  left  undisturbed  by  SURFPL0T . 
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The  plot  returned  by  SURFPL0T  is  the  same  in 
principle  as  shading  each  bounded  region  in  a 
contour  plot  with  a  distinct  shade  of  gray 
that  corresponds  to  the  levels  that  have  de¬ 
fined  the  region.  There  are  normally  ten 
distinguishable  shades  of  gray  (including 
unmarked  paper)  that  have  been  created  with 
the  line  printer  character  sets.  This  number 
is  not  ironclad  and  is  subject  to  vagaries  of 
paper  and  ribbon.  However,  the  gray  shades 
generally  are  satisfactory  using  the  full  ten 
shades  of  gray  during  daytime  runs  at  NRL  with 
either  system. 

The  operation  of  SURFPL0T  is  ordinarily  as 
follows : 

(1)  The  array  of  f(x,y)  values  is  searched 
for  the  maximum  and  minimum  values. 

(2)  The  range  from  maximum  to  minimum  is 
divided  by  ten  partitioning  values  using  a 
linear  scale  (see  Figure  1). 

(3)  The  array  of  f(x,y)  values  is  then 
searched  point-by-point  generating  a  quantized 
representation  of  each  row,  a  row  at  a  time, 
using  the  partitioning  from  step  (2). 

(4)  The  points  are  printed,  a  row  at  a  time, 
using  the  shade  of  gray  unique  to  the  quanti¬ 
zation  at  each  (x,y)  point. 

The  physical  plot  of  the  array  of  f(x,y)  values 
is  printed  so  that  the  columns  of  the  array  are 
parallel  to  the  sprocket  holes  in  the  line 
printer  paper.  The  rows  are  printed  across  the 
page.  The  line  printer  has  a  resolution  of  ten 
characters  per  inch  across  the  page  and  six 
lines  per  inch  along  the  sprocket  holes.  Hence, 
the  aspect  ratio  is  not  one.  There  is  no 
limitation  in  the  number  of  rows  (equivalent  to 
printer  lines)  other  than  the  one  provided  by 
the  line  count  on  the  RUN  card.  However,  no  more 
than  128  columns  can  be  printed  on  any  one  printer 
page.  For  arrays  that  have  more  than  128  columns, 
the  plot  is  printed  out  in  blocks  of  100 
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Fig.  1  -  Partitioning  of  Data  Range  with  a  Linear  Scale 


columns  per  page  (the  routine  automatically 
accomplishes  this  step)  with  each  page  annotated 
(see  Section  3.8). 

There  are  three  functional  pieces  to  SURFPL0T: 

(1)  generation  of  the  quantization  levels 
used  to  represent  f(x,y)  and  the  equivalent 
shade  of  gray, 

(2)  generation  of  the  annotation,  and 

(3)  generation  of  the  line  printer  image  of 
the  array  of  f(x,y)  values  with  the  notation  into 
line  printer  statements  and  the  printing. 

The  user  can  be  involved  heavily,  partially,  or 
not  at  all  in  the  first  area  and  thereby  exercise 
control  over  the  nature  of  the  plot.  The  number 
of  shades  of  gray  and  the  type  of  scale  are  iust  two 
of  the  variables  that  can  be  controlled,  if 
desired.  Only  when  the  plot  is  done  piecemeal  is 
the  user  involved  in  (3)  and  then  only  in  the  most 
minimal  fashion. 

2.2  Problem  Background 

Most  of  the  available  routines  for  plotting  in 
the  RCC  library  are  for  functions  of  one  variable. 
Representative  routines  that  can  plot  an  array  of 
f(x,y)  values  include  the  CalComp  THREE-D  Plotting 
Routine,  J5-NRL-PL0TIS0M,  J 5-NRL -SURFACE ,  and 
J5-0RNL-C0NT0UR .  These  four  routines  can  give  an 
excellent  representation  of  the  function  surface, 
however  they  require  the  use  of  a  CalComp  plotter. 
SURFPL0T  is  a  simple  plotting  routine  that  returns 
a  usable  image  of  the  function  surface  with  minimum 
cost  in  turn-around  time  and  minimal  user  involve¬ 
ment.  For  situations  where  the  point  by  point 
generation  of  an  image  of  f(x,y)  is  useful,  more 
flexible  control  of  the  image  generation  is  avail¬ 
able. 
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3 . 0  USAGE 


3 . 1  Calling  Sequence  or  Operational  Procedure 

CALL  SURFPL0T  (ZDATA,  IXLENGTH,  JYLENGTH, 

ITITLE ,  ISCALE,  ZSCALE,  ZLEVEL,  ZPARAMTR , 

IEXTEND,  ISEARCH) 

The  calling  sequence  used  is  the  same  form 
whatever  the  level  of  control  being  exercised 
by  the  user.  The  number  of  parameters  that 
need  be  specified,  however,  vary  depending  on 
(1)  a  single  or  multiple  call  to  SURFPL0T  itkI  (2)  the 
degree  of  control  exercised  by  the  user  over  the 
nature  of  the  plot.  The  user  always  names  the 
array  of  f(x,y)  values,  its  dimensions  and  sets 
up  a  title.  All  that  need  be  done  thereafter 
is  to  dimension  properly  the  remaining  arrays 
and  a  linearly  scaled  plot  will  result.  To 
exercise  control  over  the  plot,  the  parameters 
used  are  ISCALE,  ZSCALE,  and  occasionally 
ZLEVEL  and  ZPARAMTR.  To  plot  piecemeal  or  with 
several  calls  one  appropriately  sets  ZPARAMTR, 
IEXTEND,  AND  ISEARCH. 

3 . 2  Arguments,  Parameters,  and/or  Initial  Conditions 

The  user  must  dimension  all  arrays  in  the  calling 
program. 

ZDATA  -  the  name  of  the  array  containing  the 
f(x,y)  values,  which  can  be  either 
fixed  or  floating  point.  The  array 
is  left  undisturbed  by  SURFPL0T. 

DIMENSION  ZDATA  (IXLENGTH,  JYLENGTH) 
IXLENGTH  -  the  number  of  rows  in  ZDATA 
JYLENGTH  -  the  number  of  columns  in  ZDATA 
ITITLE  -  the  name  of  the  array  containing  the 

plot  title,  one  card  of  BCD  characters. 
DIMENSION  ITITLE  (10).  See  Section  7.0. 

While  the  above  four  parameters  must  be  specified 
in  each  call,  the  remaining  parameters  are  optional. 
The  next  group:  ISCALE,  ZSCALE,  and  occasionally 
ZLEVEL  and/or  ZPARAMTR,  are  used  to  control  the 
plotting . 
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I SCALE 


-  the  name  of  an  array  of  four  numbers 
that  selects  the  number  of  shades  of 
gray,  the  type  of  scale,  a  "positive" 
or  "negative",  and  indicates  whether 
the  elements  of  EDATA  are  fixed  or 
floating  point  numbers.  DIMENSION 
ISCALE (4) . 


0 


ISCALE(l)  = 


ten  shades  of  gray 
(as  shown  in  the  sample 
plot  in  Section  7.0.) 


1-25  ISCALE(l)  shades  of 
gray  -  if  more  than 
ten,  the  user  must 
supply  the  symbol  set 
using  three  data  cards. 
These  cards  are  read  in 
under  a  25(R1)  format. 
There  is  an  interaction 
with  ISCALE(3) ;  when 
zero  the  first  symbol 
on  the  card  will  repre¬ 
sent  the  smallest 
amplitude  and  when  non¬ 
zero  the  largest. 

See  Section  3.7. 


ISCALE ( 2) 


3 

\4 


Linear  Scale 
Logarithmic  Scale 
Logarithm  of  sine  to 
the  pth*  power 
sine  to  the  pth*  power 
user  supplied;  the 
quantizing  levels  must 
be  in  ZLEVEL 
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increasing  amplitude 
printed  as  increasing 
darkness  -  a  "negative". 
Probably  the  most  useable 
form. 

decreasing  amplitude 
printed  as  increasing 
darkness  -  a  "positive". 


*If  ISCALE (2)  =  2,  p  =  Z.5 
*If  ISCALE (2)  =  3,  p  =  0.25 


7 


ISCALE  (4) 


0  -  the  elements  of  the  array 
ZDATA  are  floating  point 
numbers . 

1  -  the  elements  of  the  array 
ZDATA  are  fixed  point 
numbers . 


ZSCALE  -  the  name  of  an  array  of  three  numbers 

used  to  generate  the  quantizing  levels  in 
ZLEVEL.  A  full  explanation  of  the  meanings 
of  ZSCALE  is  deferred  to  Section  4.0  where 
the  algorithms  for  generating  ZLEVEL  are 
explained.  Along  with  ISCALE,  ZSCALE  can 
be  used  to  exercise  complete  control  over 
the  plotting.  By  leaving  ZSCALE(*)  =  0.00, 
the  present  parameters  will  be  selected. 

To  use  ZSCALE(*)  with  other  than  nominal 
values,  see  Tables  1  §  2,  Section  4. 
DIMENSION  ZSCALE (3). 


ZLEVEL  -  the  name  of  an  array  of  quantizing  values 
generated  by  SURFPL0T .  If  ISCALE(2)  ^4, 
the  user  need  only  dimension  ZLEVEL.  If 
ISCALE(2)  =  4,  the  user  supplies  the 
quantizing  levels  in  ZLEVEL  (see  Section  4.) 
DIMENSION  ZLEVEL (10)  or  ZLEVEL(ISCALE(1) )  if 
more  than  ten  shades  of  gray  are  used. 

The  last  group  of  parameters  is  used  when  the  plot  is 
made  with  several  successive  calls  to  SURFPL0T. 


ZPARAMTR  -  the  name  of  an  array  containing: 


ZPARAMTR(l) 
ZPARAMTR (2) 

ZPARAMTR  (3) 

ZPARAMTR  (4) 
ZPARAMTR  (5) 

ZPARAMTR  (6) 


=  maximum  value  of  f(x,y) 

*  value  of  x  for  f  (x,v) 

-  value  of  y  for  f  (x,v) 

=  minimum  value  of  f(x,y) 

=  value  of  x  for  f  .  (x,y) 

=  value  of  y  for  f  .  (x,y) 

1  minv 


In  the  multiple  call  situation,  to  have 
the  same  scale  for  all  portions  of  the 
plot,  the  scale  routines  must  be  given 
the  maximum  and  minimum  values  of 
f(x,y)  at  the  start.  This  is  accomplished 
by  filling  ZPARAMTR(')  with  the  appropriate 
values.  If  left  zero,  the  first  portion 
of  the  array  that  is  plotted  will  be  search¬ 
ed  and  the  scales  from  this  first  portion 
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will  be  used  for  the  remaining 
portions.  In  the  single  call  case, 
ZPARAMTR  returns  the  maximum  and 
minimum  parameters  of  the  plotted 
array.  DIMENSION  ZPARAMTR (6) . 

I EXTEND  -  the  name  of  an  array  containing  multiple 
call  indicators;  IEXTEND(l)  denotes  that 
the  extensions  will  be  in  the  x-direction 
(parallel  to  the  sprocket  holes),  IEXTEND(2) 
for  y-direction  extensions.  Either  one  or 
both  extensions  can  be  called.  DIMENSION 
IEXTEND(2) .  The  following  code  is  used. 
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IEXTEND  = 


only  one  call  in  this 
direction 

first  of  successive  calls 
in  this  direction 
neither  first  nor  last 
call 

last  of  successive  calls 


ISEARCH  -  a  variable  used  in  multiple  call  situations 
to  force  the  scales  to  be  recomputed  for 
each  separate  call  in  the  sequence.  This 
means  that  the  scales  change  (presumably) 
for  each  call.  The  only  reason  for 
including  this  option  is  that  the  maximum 
and  minimum  values  of  f(x,y)  can  be  found 
for  the  extended  array  from  the  maximum 
and  minimum  values  returned  in  ZPARAMTR(-) 
for  each  section  of  f(x,y).  A  p^t  of 
each  section  also  is  printed. 


ISEARCH 


!0  -  scales  are  not  recomputed 
1  -  scales  are  recomputed 


3.3  Space  Required  (Decimal  and  Octal) 

3.3.1  Unique  Storage 

6100  octal  (3136  decimal)  locations 

3.3.2  Common  Blocks 
None 

3.3.3  Temporary  Storage 
None 
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3 . 4  Messages  and  Instructions  to  the  Operator 
None 

3 . 5  Error  Returns,  Messages,  and  Codes 
None 

3 . 6  Informative  Messages  to  the  User 
None 

3 . 7  Input 

In  addition  to  a  card  containing  the  title,  input 
is  required  when  the  user  supplies  the  symbol  set 
to  define  the  shades  of  gray.  Up  to  25  shades  can 
be  defined.  Each  shade  is  made  up  of  three 
characters  from  the  line  printer  set  as  read  in 
from  three  data  cards  punched  in  a  2 5 (R1 )  format. 
In  the  example  below  the  first  gray  level  is  the 
character  A  overprinted  with  a  B  and  C. 


3 . 8  Output 

A  typical  output  is  shown  in  Section  7.  In  addition 
to  the  line  printer  plot,  the  following  information 
is  printed: 

(N,M)  SURFACE  PLOT . Page  N1  of  N2 

The  (N,M)  refers  to  multiple  call  situations 
and  is  used  to  annotate  each  portion  of  the 
overall  plot  as  an  element  in  a  matrix.  N  is 
the  row  i:\dex  and  M  is  the  column  index.  For 
the  single  call  case,  N  =  M  *  1 . 

The  page  numbering  accounts  for  arrays  with 
more  than  128  columns;  then  the  plot  will 
be  made  100  columns  per  page  with  the  N1  and 
N2  annotating  this  string. 
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ITITLE  -  the  user  supplied  title. 


DATA  MAXIMUM . the  maximum  value  of  f(x,y)  and 

its  (x,y)  location  in  the  array 
(the  first  occurrence  of  ^max » y ) ) 

DATA  MINIMUM . the  minimum  value  of  f(x,y)  and 

its  (x,y)  location  in  the  array 
(the  first  occurrence  of  fm^n(x,y)) 

DATA  RANGE . if  linear  plot,  (fmax-fmin)  J  if 

logarithmic,  C-log10(fmax/fmln) 


where  C  is  equal  to  10.00  unless 
otherwise  specified  by  the  user. 
See  Tables  1  §  2,  Section  4. 


SCALE  MAXIMUM . maximum  value  of  the  quantizing 

levels 

SCALE  MINIMUM . minimum  value  of  the  quantizing 

levels 

SCALE  RANGE . if  linear,  (Scalemov  -  Scale  .  ) 

’  v  max  min 

quantizing  levels;  if  log, 

C'logio  (Scalemax/Scalemin’- 

C  defined  as  above. 

PLOT  REFERENCE . if  linear,  (f  -f  J  unless 

max  nun 


specified  by  the  user;  if  log, 

either  f  ,  f  /2.00,  or  user 
max ’  max  ’ 

supplied.  See  Tables  1  f,  2 , 

Section  4.0. 


PLOT  BIAS . for  linear  scales,  nominally 

f  .  (x,y)  unless  otherwise 
minv  ,7J 

specified  by  the  user;  if  log 
scale,  it  will  be  0.00. 

PLOT  EXTENT . the  dimensions  of  the  array  of 

f(x,y)  values  ( I XLENGTH ,  JYLENGTU) 

DB  REFERENCE . the  "C"  used  to  compute  decibels. 

If  ZSCALE (2)  =  0,  DB  REFERENCE  = 
10.00;  if  ZSCALE (2) /0,  DB  REFERENCE 
ZSCALE  (2)  . 
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HALF  TONE  DENSITY  WEDGE  AND  SCALE 


The  equivalence  between  shades  of 
gray  and  numerical  values  and  the 
relation  to  the  data.  If  linear, 
percentage  of  the  maximum;  if  log, 
decibels  relative  to  the  maximum*. 
Five  rows  of  the  gray  scale,  each 
gray  shade  repeated  for  ten  points, 
are  printed.  If  more  than  ten 
shades  of  gray  have  been  specified 
the  wedge  is  printed  along  the 
page  rather  than  across  the  page. 

The  body  of  the  plot  follows.  A  border  of  M  is 
printed  around  the  periphery.  Every  fifth  row  and 
column  are  numbered,  the  maximum  and  minimum  of  the 
array  of  f(x,y)  values  are  labeled  (only  the  first 
occurrence  of  each).  The  program  will  keep  track  of 
multiple  call  plots  so  that  the  numbering  is 
consistent  and  the  borders  removed  appropriately  to 
ease  the  cutting  and  pasting. 

After  completion  of  the  call,  the  array  ZPARAMTR 
contains  the  maximum  and  minimum  parameters  of  the 
array  (see  Section  3.2)  and  the  array  ZLEVEL 
contains  the  quantizing  level  values. 

3.9  Formats 


See  Section  3.7. 

3.10  External  Routines  and  Symbols 


None 

3.11  T  iming 

The  routine  takes  roughly  0.07  seconds  per  row  of 
the  array  of  f(x,y)  values. 


*decibels  are  conventionally  defined  as  10  log(-)  =  dB 


3.12  Accuracy 

Not  applicable 

3.13  Cautions  to  Users 

The  aspect  ratio  of  the  plot  is  not  one  so  a 
distortion  exists  due  to  the  stretching  of  the 
x-axis  relative  to  the  y-axis. 

This  imaging  does  not  work  well  with  small  or 
thin  arrays,  for  the  technique  depends  on  the 
eye  perceiving  structures  in  the  data.  If  the 
gray  shade  does  not  have  enough  spatial  extent, 
the  eye  does  not  perceive  a  pattern. 

3.14  Program  Deck  Structure 

gJ0B 

Iftn 

user  program  (contains  call  to  SURFPL0T) 
SUBROUTINE;  SURFPL0T 
SC0PE 

gL0AD 

gRUN 

data  card  1  -  title  card 

data  card  2|  plotting  symbols  if  more  than 
data  card  3>-  ten  shades  of  gray  are 
data  card  4;  requested, 
end-of-file  card 

3.15  References  -  Literature  -  Appendices 

[1]  D.  E.  Knuth,  The  Art  of  Computer  Programming, 
Volume  2  -  Seminumer ical  Algorithms, 

Addison  -  Wesley,  1969. 
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Fig.  2  -  Rough  Transfer  Characteristics  for  Data  Scales  0  through  3 


4.0  METHOD  OR  ALGORITHM 


Representing  a  function  of  two  variables  on  a  flat 
piece  of  paper  so  that  its  structure  is  obvious  to 
the  eye  requires  some  trickery.  The  surface  can  be 
suggested  by  encoding  the  amplitude  of  each  (x,y) 
point  into  a  pseudo-height  interprctable  by  the  eye. 

Some  standard  methods  are: 

(a)  mapping  amplitude  into  brightness  as  in  a 
photograph , 

(b)  mapping  constant  amplitude  contours  as  in 
the  topographic  map, 

(c)  using  isometric  distortions  (or  projections) 
by  mapping  amplitude  for  each  row  (or  column) 
into  y-axis  deflection  then  offsetting  each 
succeeding  row  a  fixed  amount  in  the  x-and 
y-d irect ions ,  and 

(d)  mapping  regions  of  constant  amplitude  or 
amplitude  bands  into  colors  or  shades  of 
gray  as  in  a  color  relief  map. 

SURFPL0T  follows  method  (d) .  The  data  f(x,y)  is 

partitioned  into  bands,  each  band  is  assigned  a 

unique  shade  of  gray,  and  the  array  of  f(x,y)  values 
is  printed  point  by  point  following  this  code. 

The  technique  of  selecting  the  amplitude  bands  is  not 
straightforward.  In  Figure  2  four  different  parti¬ 
tionings  of  the  data  range  are  displayed.  The 
attraction  of  non-linear  scales  is  simply  that  it 
allows  the  resolution  of  the  eye,  represented  by  the 

ten  levels  along  the  abscissa,  to  be  placed  over  the 

significant  detail  of  the  data.  A  classical  example 
is  a  monochromatic  image;  the  array  of  f(x,y)  values 
representing  the  image  often  has  considerable  dynamic 
range.  However,  on  closer  examination,  the  distribution 
of  amplitudes  often  turns  out  to  be  confined  to  a  small 
portion  of  the  dynamic  range  not  far  from  the  minimum 
value  with  a  few  very  large  specular  returns.  The 
information  content  is  in  this  small  portion  of  the 
dynamic  range  near  the  minimum;  hence  a  logarithmic 
scale  would  be  an  appropriate  technique  to  use  for 
partitioning  the  dynamic  range,  (see  Figure  2). 
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In  constructing  a  partition  of  the  data  there  are 
three  features  available  to  the  user  for  control: 

(1)  the  number  of  bands, 

(2)  the  shape  of  the  partitioning  curve,  and 

(3)  the  location  of  the  partitioning  curve. 

Examples  of  (2)  are  shown  in  Figure  2.  An  example  of 
(3)  for  linear  scales  is  shown  in  Figure  3.  In 
SURFPL0T,  the  first  two  features  are  specified  by 
ISCALE  (1)  and  ISCALE  (2)  respectively.  To  specify 
the  third  feature,  the  program  uses  ISCALE(l), 
ZSCALE(l),  ZSCALE(2)  and  ZSCALE(3).  The  partitioning 
values  are  located  in  array  ZLEVEL.  Table  1  shows  how 
each  scale  is  computed;  Table  2  shows  the  specification 
of  the  third  feature.  The  following  pattern  is  used; 
if  ZSCALE(-)  is  left  zero,  a  preset  value  is  used.  If 
non-zero,  that  value  of  ZSCALE(*)  is  used.  In  the  case 
of  ZSCALE(l),  a  third  option  is  provided.  If 
ZSCALE(l)  =  -1.00,  then  half  the  preset  value  is  used. 
This  is  most  useful  in  logarithmic  plots  as  it  slides 
the  scale  down  3dB  below  the  maximum. 

The  plotting  consists  of  printing  one  of  a  series  of 
half- tone-like  shades  of  gray  at  each  point  correspond¬ 
ing  to  the  f(x,y)  value  at  that  point.  The  shade  of 
gray  is  determined  by  finding  the  partitioning  level 
that  just  equals  or  exceeds  the  f(x,y)  value.  Assigned 
to  each  partitioning  level  is  a  unique  shade  of  gray 
created  by  overprinting  wit i  three  line  printer 
characters.  The  array  is  examined  row  by  row;  printing 
is  also  accomplished  row  by  row.  The  data  array  is 
unaltered. 
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TABLE  1  -  SCALE  ALGORITHMS 
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TABLE  2  -  DEFINITION  OF  TERMS 


40 
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5.0  SOURCE  LANGUAGE  LISTING 


JS 


C 

C 

C 

C 

C 

C 

C 

C 

C 

r 

c 

c 


NRL  SURFPI  OT  LINF  PRINTER  P|01S  OF  TW0-DTMFNS10NAL  SURFACFS 

SUBROUTINE  SIIRFPL  OT  (70  at  a.  f  XLFNjfiTH.  JYLFNGTH,  I T I TLF  » ISCALF»2SCALE  * 
1  7LFVFL ,7P A PAMTR, l FX TEND, TSF ARCH) 


IOFNT  NIIMRFW  -  J5008000 

TITIF  -  L INF  PRINTER  PLOTS  OF  TWO  DIMENSIONAL  SURFACFS 
IOFNT  NAMF  -  J5-NRI -SIIRFPLOT 
I  ANGtJAGF  -  FORTRAN 

computfr  -  cnr-moo 

CONTRIBUTOR  -  MICHAFL  A.  TAMNY,  COnF  5365  T «  A I RROPNF 

RAOAR  DIVISION 


nooooioo 
00000200 
00000300 
00000400 

ooooosoo 
00000600 
00000700 
00000800 
RADAR  BRANCH, 00000900 
00001000 


ORGANISATION  -  NR| 


NAVA)  RFSFARCH  LARORATORY  - 


WASHINGTON, 

20375 


D.C. 


OATF  -  T I  OCTOHFR  1973 

PtJRPOSF  -  TO  PLOT  AN  ARRAY  OF  F(X.Y)  VALUFS  USING 
RY  RFPRFSFNTING  AMPl I TUDF  WITH  A  SHADF  OF  GRAY 


DIMENSION 
niMFNSION 
DIMFNSION 
D  I MENS  I  ON 
OIMFNSION 
OIMFNSION 
DIMFNSTON 
DIMFNSION 
OIMFNSION 
DIMENSION 
0 IMFNS  T  ON 


70ATA  I  1  ) 

7PARAMTR  (  1 ) 

IFXTFNOI  1 > 

ISFAI  Fill 
7SCALF  I  1  ) 

7LFVFL  (  1  ) 

ititi  fm  > 

I  1NFPI  OT(13S,3) 

I  INF1  I  1 3S ) ,  )  INF? 1 1 3S ) , 
.11  INF  117,3) 

)  I  NFS YM I ?S , 3 ) 


I INF31135) 


00001100 
00001?00 
00001300 

THE  LINF  PRINTFR00001400 
00001500 
00001600 
00001700 
00001800 
00001900 
00002000 
00002100 
0000??00 
0000?300 
00002400 
00002500 
00002600 
00002700 
00002800 
00002900 
00003000 
00003100 
00003200 
00003300 
00003400 
00003500 
00003600 
00003700 
00003800 
00003900 
00004000 
00004100 
00004200 
00004300 
00004400 
00004500 

•,1R*,1R*,1RX,?I1R0)  ,IR|),1R0,1RD)  00004600 
DAT  A  (  II  TNFSYMI !.»?)«!  =I«10)=6(1P  I  ,  1RX,  1PA  ,  1R»  •  1RP)  00004700 

DATAIO  TNFSYMIl ,3),l =l»10)s8(lP  >,1RU,1RX)  00004800 

DATA  (TCHIFTl  =  1 OOOOOOOOOOOOnOR)  00004900 


FOIIIVAI  FNCF 
FQUTVAI  FNFF 
FOIIIVA)  FNCF 
FOII I  V«|  FNCF 
FOIIIVAI  FNCF 


I.INFPLOTI  1 )  ,L  INF  1(1) 
I  INFPLOT  1 1  36)  ,|  INF? II) 
I  1NFPI  OT (271 ) »L  INF31 1 ) 
I  INFPLOT  (  1  )  ,  J|  INF(l) 

7TFST , I  TEST  ) 


TYPF 

TYPF 

TYPF 


intfgfr 

intfgfr 

INTFGFR 


0)  NKF1I 
Rl  ANK. 
AF  II  I.  , 


L,  ORF ILL 
DOT,  P|  US 
IFILI  ,  MFILI 


,  NFILI  ,  UFU  I  ,  WF  ILL  ,  XFILL 


rata  (RLNKF  ILL  =  ?R  ),  lOMFIIL 


DATA 

( R|  ANK 

r 

1R  )  , 

(OOT  = 

1R 

• ) , 

DATA 

(AFII  1 

= 

IRA)  , 

(IFILI 

= 

1  R I ) 

OATA 

(WFILI 

= 

1 PW )  , 

(XFILL 

= 

IPX) 

DATA 

(NX  = 

I  )  i 

,  (NY  = 

1  1 

DATA  ( 

(1  INFSYMi 

ll,l),l 

=1,10)= 

:  1 R 

,  1 R 

2ROR) 

(PLUS  =  1 R  ♦  ) 

,  IMF  ILL  *  1RM). 


INF  II  L  =  IRN) 
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n  r>  r>  o  o  o  r>  ooooo 


4  / 


DATA 

OATA 

DATA 

OATA 

OATA 

OATA 


(ISHIFT2 
(TSHIFT3 
(ISHIFT4 
( ! SHIFTS 
(  T5HIFT6 
< ISHIFT7 


1000000000000R) 

10000000000B) 

1 00000000B) 

1 OOOOOOB) 
ioooor) 

1000) 


GO  TO  (  1005*1000  )  NX  ♦  NY 


THIS  TRAP  15  USED  TO  RECOGNIZE  THE  T YPF  CF  CYCLF  IN  EXTFNOED 
CALL 5  TN  ROTH  DIMENSIONS.  THF  FIRST  CALL  OF  THIS  TYPF  CAUSFS 
NX  =  NY  =  0.  OTHERWISE  THEY  ARE  INOICIFS  USFD  TO  KFEP  TRACK  OF 
THF  ARRAY  BLOCKS  SENT  TO  SURFPLOT  FOR  PRINTING  AND  AS  SUCH  ARE 
GREATER  THAN  OR  FOUAl  TO  1. 

100R  IF  (  IFXTFND<1)-IEXTFN0<?)  .EO.  1  )  GO  TO  1006 
NX  =  1 
NY  -  2 
ILAST  =  0 
IOXY  =  2 
GO  TO  1000 

1006  NY  *  1 
NX  =  2 
-ILAST  =  0 
IOXY  -  1 

THF  ARRAY  IS  SFARCHFU  FOR  THF  MAXIMUM  AMD  MINIMUM  VALUES  AND 
THF  SCALFS  COMPUTED  IN  THE  FOLLOWING  MANNER .  •  ( 1  )  IN  SINGLE  CALL 
CASES.. I?)  IF  I SE ARCH  =  1..H)  TF  THE  FIRST  CAIL  OF  AN  FXTENOEO 
CAI  L  ANn  THE  MAXIMUM  ANO  MINIMUM  VALUFS  HAVE  NOT  BEEN  PASSED  IN 
7PARAMTR.  IN  THE  SF ARCH*  LOCATIONS  ARF  NUMBERFO  IN  THE  FXTENOEO 
ARRAY  BY  BIASING  WITH  ILAST  ANO  JL  AST •  IF  A  SF ARCH  IS  MADE  THF 
LEVELS  MUST  RF  COMPUTFO. 

100«  I F (  IFXTFNOIl ) *IFXTFN0<2)  ,F0.  0  )  GO  TO  1009 
IF (  ISFARCH  ,F0.  1  )  GO  TO  1009 

IF (  IFXTFNO (1)  .GT.  1  .OR.  1FXTFN0I?)  .GT.  1  )  GO  TO  1075 

IF (  ZPARAMTR(l)  .FO.  0.00  .AND.  7PARAMTR (4)  .FO.  0.00  )  GO  TO  1 

GO  TO  10?0 


1009  7TFST  =  ZOATAll  ) 

IF (  ISCA|F(4)  .FQ.  0  )  GO  TO  1010 
7TFST  =  ITFST 

1010  7  MAX  =  7  MIN  =  7TFST 
MAX  X  =  MIN  X  =  l  ♦  ILAST 


00005000 
00005100 
00005200 
00005300 
00005400 
00005500 
00005600 
00005700 
00005A00 
00005900 
00006000 
00006100 
00006110 
00006120 
00006130 
00006140 
00006150 
00006200 
00006300 
00006400 
00006500 
00006600 
00006700 
00006800 
00006900 
00007000 
00007100 
00007200 
00007300 
00007400 
00007410 
00007420 
00007430 
00007440 
00007450 
00007460 
00007470 
00007500 
00007600 
00007700 
00007800 
00900007900 
00008000 
00008100 
00008200 
00008300 
00008400 
00008500 
00008600 
00008700 
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4  ► 

4  * 


n  n 


MAX  Y  =  M  TN  Y  =  1  ♦  JL  AST 
00  IOIS  J  =  1 • JYLENGTH 
INOEX  =  (  J-l  I  *  I  XLENGTH 
no  IOIS  T  =  1*1 XLENGTH 
7TFST  =  ZOATA ( TNREX ♦  I  > 

IF  (  ISCALF (41  .FO.  0  )  GO  TO  1011 
7TEST  =  ITFST 

1011  IF  (  7TFST  .LF.  7  MAX  I  GO  TO  101? 

7  MAX  =  7TFST 
MAX  X  =  I  ♦  HAST 
MAX  Y  =  J  ♦  JL AST 
GO  TO  IOIS 

101?  IF  (  7TFST  .GF.  7  MIn  >  GO  TO  IOIS 
7  MIN  =  7TFST 
MIN  X  =  I  ♦  HAST 
MIN  Y  =  J  ♦  JL AST 
IOIS  CONTINUE 

7PARAMTR ( 1 )  =  7  MAX 
7P ARAMTR ( ? 1  =  MAX  X 
7PARAMTR ( 1)  =  MAX  Y 
7P AR AMTR (4 )  =  7  MIN 
7PARAMTR (S)  =  MTN  X 
7P AH AMTR  ( 6 )  r  MIN  Y 

IF  (  IFXTENO ( 1 )  .GT.  1  .OR.  !FXTENn<?»  ,GT.  1  )  GO  TO  10S0 
GO  TO  1030 

10?0  7  max  =  7PARAMTRM) 

MAI  X  =  7PARAMTR (?) 

MAX  Y  =  7pAHAMTR ( 3) 

7  MIN  =  7PARAMTRI4) 

MIN  X  =  7P ARAMTR ( S ) 

MIN  Y  =  7P ARAMTR ( 6 ) 

this  section  sets  the  number  oe  gray  shares.  and  if  less  than 
TFN»  SFLFCTS  ones  FURTHER  APART  USING  I  STEP  AmO  LRTAS. 

1030  IE  (  1  SCALE ( 1 )  .GT .  10  »  GO  TO  1040 

IE  (  I  SCALE ( 1 )  .FO.  0  )  IlF/ELS  =  10 

IF  (  ISCALE(l)  .NF.  0  )  I  LEVFl S  =  ISCAlE(I) 

LSTFP  r  10/LLFVELS 

LB  IAS  =  10  -  LSTEP*LLEVELS 

GO  TO  JOSO 

1040  LLEVFLS  =  1SCAI  Ed) 

LSTFP  =  1 
L  B I A  S  =  0 

REAR  91041.  (  L  INFSYM  (I  .  I  )  .  L  =  1  .11  EVFI.S  ) 

REAR  Q1  04 1  »  (I  INFSYM  (I  ,?)  ,  L  =  1  .1  I  FVFl.S  ) 

REAR  01041.  (  I  INESYMU  ,1),  L  -  I.ILEVFLS  ) 

<51041  FORMAT!  ?S  (  R1  )  ) 

C  THIS  SECTION  COMPUTES  THE  PARTITIONING  IFvFLS.  7LEVEL 


OOOORROO 
00008900 
00009000 
00009100 
00009?00 
00009300 
00009400 
00009S00 
00009600 
00009700 
00009R00 
00009900 
00010000 
00010100 
00010?00 
00010300 
00010400 
00010500 
00010600 
000)0700 
0001 OROO 
00010900 
00011000 
00011100 
00011200 
00011300 
00011400 
OOOllSOO 
00011600 
00011700 
0001 1R00 
00011900 
00011910 
00011920 
00012000 
00012100 
00012200 
00012300 
00012400 
00012500 
00012600 
00012700 
0001 ?R00 
00012900 
00013000 
00013100 
00013200 
00013300 
00013400 
00013410 


1050  GO  TO  (  10S1 ♦10SS»I060,|065.|070  )  TSCAI.FI?)  ♦  1 


10S1 


ITYPF  =  RLNKFII.L 
IF(  /SrULFIl)  .FO.  0.00 
IF (  7SFA1  F<1>  .F0.-l.00 

IF (  7SFALFM)  .NF.  0.00 
7  RFF  =  ZSFALFll) 

IF  (  7SFALFI?)  .FO.  0.00  )  / 
IF  (  7SFALF<?)  . NF .  0.00  )  7 
ORSFAIF  =  0.00 
STFP  =  7RFF/ILFVFLS 
00  105?  L  =  1  .11  FVFLS 
7LFVFLH  »  =  7  HIAS  ♦  STFP*| 
FONT  l Nl IF 

7RANGF  r  7  MAX  -  7  MIN 
SCALF  =  7LFVFL (LLFVFI  S> 

GO  TO  1075 


10SS  ITYPF  =  OOF  II.  I. 


> 


)  7  RFF  =  (  7MAX-7MIN 
)  7  RFF  =  (  7MAX-7M1N  >/?.00 
.ANO.  7SCALFin  .NF.-l. 00  > 


0  T  AS  = 
0  I  AS  = 


7  MIN 
7SCAI.F1?> 


10S? 


-  7LFVFI  (1> 


10SG 


STFP  =  1 .00/LLFVFl  S 

no  loss  I.  =  1  .LI  FVFLS 

POWFR  =  -SI OPF * (  I .  00-STFP*L  > 

7LFVFLU)  =  7RFF  *  (  10.00**POWFR  > 

FONT  IN'IF 

7RANGF  =  DHSFAl  F»Al  OGI 0  (  7MAT//MIN 


) 


SCA|  F  =  nosCA|  F*A|  OGI  0  ( 
GO  TO  1 0 7S 

1060  1TYPF  =  OOF  ILL 


7\  FVF|  ULFVFlS>/7LFVFMI»  > 


ooonsoo 
00011600 
00013700 
00013000 
00013900 
00014000 
00014100 
000 14?00 
00014300 
00014400 
00014S00 
00014600 
00014700 
00014R00 
00014900 
00015000 
0001S100 
000  I 5?00 
000 1 5300 
00015400 


IF  < 

7SFALF < 1 > 

.FO. 

0.00 

) 

7  RFF  = 

7  MAX 

00015S00 

IF  ( 

7SFALF ( 1 > 

.FO. 

-1.00 

> 

7  RFF  ^ 

7MAX/P.00 

00015600 

IF  ( 

7SF  Al.F  (  1  > 

.NF. 

0.00 

• 

ANn.  7SFAI  F ( 1 )  . NF.-l .00  ) 

0001S700 

7  RPF  =  7SFA| F  <  1  ) 

00015R00 

IF  ( 

7SF  ALF (?) 

.FO. 

0.00 

> 

nRSFAl  F 

=  10.00 

00015900 

IF  ( 

7SFAI  F(?l 

.NF. 

0.00 

> 

nRSFALF 

=  7SCALFC?) 

00016000 

IF  ( 

7SFALF 11) 

.FO. 

0.00 

) 

si  opf  = 

40 . OO/DPSCAI  F 

00016100 

IF  ( 

7SFALF ( 3) 

.NF. 

0.00 

> 

SLOPF  = 

7SCAI  F  (  3) 

0001 6?00 

7  H  f  AS  =  0.00 

00016300 

00016400 
00016500 
00016600 
00016700 
00016000 
00016900 
00017000 
00017100 
0001 7?00 
00017300 


IF  ( 

7SF At  F ( 1 )  .FO.  0.00 

> 

7  RFF  = 

7  MAX 

0001  /400 

IF  ( 

7SF  ALF ( 1 )  .FO.-l.nO 

> 

7  RFF  = 

7MAX/P.00 

00017500 

IF  ( 

7SFA|  F ( 1 )  .NF.  0.00 

• 

ANO.  7SFAI  F 11  >  .NF.-l .00  > 

00017600 

7 

RFF  =  7SFAI  F  11) 

00017700 

IF  ( 

7SFAIF(?>  .FO.  0.00 

> 

nBSFAl  F 

=  10.00 

00017R00 

IF  < 

7SF  ALF (?)  .NF.  0.00 

> 

OHSFAI  F 

=  7SFALF (?) 

00017900 

IF  I 

7SFAI  F  (  3 )  .FO.  0.00 

> 

POWFR  = 

?.50 

0001R000 

IF! 

7SFAIFI3)  .NF.  0.00 

) 

POWFR  = 

7SCAI  F  13) 

000 1 R 1 00 

7  R 1  AS  =  0.00 

0001fl?00 

STFP 

=  1.00/1  t.FVFI.S 

000l?300 

PinvFW?  =  3. l4|S9PSS4/?.00 


000IF.400 


> 

! 

! 


RADIAN  =  P10VFP2«STFP 
no  1061  L  =  l.LLEVELS 
APGUMFNT  =  RADIAN*!. 

ZLFVFL  (L )  =  Z«FF*(  (SINF(ARGUMENT) >**PnwFR  ) 

1061  FONT  INlIF 

SCALF  =  nRSCALF*AL0G10{  7LFVEL(LLEVFLS»/ZLEVFL  m  ) 

7RANGF  =  ORSC  ALF*A|.OG  1  0  (  7MAX/ZM1N  ) 

GO  TO  1 n75 

1 06*5  ITTPF  *  RLNKF1LL 

IF  (  ZSCALFM)  .FO.  0.00  I  7  RFF  =  (  7MAX-7MIN  ) 

IF  (  ZSCALEO)  .F0.-l.00  )  7  RFF  =  <  ZMAX-7MIN  1/2.00 

IF  (  7SCALF  (  1  )  .NF.  0.00  .AMD.  7SCALFU)  .MF.-I.00  ) 

1  7  RFF  *  7SCALF ( I ) 

IF  (  7SCALF(?)  .FO.  0.00  )  Z  fi  I  AS  =  7  MIN 

IF  {  7SCtLF<2)  .NF.  0.00  I  7  RIAS  =  7SCALE<2> 

I F  <  7SCALF<3)  .FO.  0.00  )  POWFR  a  0.2S 
IF  (  7SCALF<3>  .NF.  0.00  I  POWFR  *  7SCAI.FO) 

DRSCALF  a  0.00 
STFP  =  1.00/LLFVELS 
PIOVFR?  a  3.141592054/2.00 
PARIAN  a  PI0VFP2*STFP 
no  1000  1  =  l.LLFVFLS 
ARGUMENT  =  RADI AN*L 

7LFVFLII)  =  7  H  I  AS  ♦  7RFFM  <S INF < ARGUMENT >> **POWFR  ) 
inoo  roNTiN'iF 

SCALF  a  7LFVFL(LLFVFLS)  -  Z1FVFL«1) 

7RANGF  =  7  MAX  -  Z  MIN 

GO  TO  1075 

1070  ITTPF  a  R|  NKFII  L 

7RANGF  =  7  MAX  -  7  MIN 

SCAlF  a  7LFVFL(Ll FVFLS)  -  ZLFVFL ( 1 1 

GRSCAI  F  a  0.00 

7  RIAS  a  7  MIN 

7  RFF  =  7  MAX 

C  THIS  SF0T10N  IS  TFST1NG  THF  ARRAY  WIDTH.  THF  LOGIC  15  TO  PRINT 

C  OUT  AS  many  OF  THF  COLUMNS  AS  WILL  FIT  ON  A  COMPUTER  PAGF  AT  ONF 

C  TIMF.  up  TO  1?«  COLUMNS  CAN  PF  PRINTED  ON  ONE  PAGF*  IF  MORE  THAN 

C  1 28  COl  UMNS  IN  ThF  ARRAY,  THEY  ARF  PRINTED  IN  PLOCKS  OF  100. 

C  An  THF  ROWS  ARF  PR  IN  TFD  IN  EACH  HLOCK. 

li>  75  JRJAS  =  Jl  AST 

IF  I  JY|  FNGTM  .GT.  1 2fl  )  GO  TO  10R0 

NPAGFS  a  1 

JLOMG  a  J Y I  tNGTH 

IF  (  JYI  FNGTH  .I  T.  114  )  JSHIFT  a  10 
IF (  J YL  FNGTH  .GT.  1)1  )  JSHIFT  a  0 

go  to  ?non 


00010500 

00018600 

00018700 

00018800 

00018900 

00019000 

00019)00 

00019200 

00019300 

00019400 

00019500 

00019600 

00019700 

00019800 

00019900 

00020000 

00020100 

00020200 

00020300 

00020400 

00020500 

00020600 

00020700 

00020800 

00020900 

00021000 

00021100 

00021200 

00021300 

00021400 

00021500 

00021600 

00021700 

00021800 

00021900 

00022000 

00022100 

00022110 

00022120 

00022130 

00022140 

00022150 

00022200 

00022300 

00022400 

00022500 

00022600 

00022700 

00022800 

00022900 
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1  ORO  NPAGFS  r  JYIFNGTH/100 

I F<  JYt  PNGTH  .GT.  NPAGFS*IOO  )  NPAGFS  =  NPAGFS  ♦  1 
JLONG  =  ion 
JSHIFT  =  10 


pooo  no  ?snn  npagf  =  I. npagfs 

C  THIS  IS  THF  PP  IN  f  I NG  SECTION 

c 

C  THIS  SECTION  SFLFCTS  THF  PROPER  SLOCK  OF  COl  UMNS  FOR  PRINTING? 

C  POSITIONS  THF  PLOTTING  ON  THF  PAGF.  ANO  SHI  CTS  THF  TYPF  OF 

C  ANNOTATION.  A  HFAOFP  IS  At  WAYS  PRINTEO  FOR  THF  FIRST  CALL  AND 

C  FXTFNDFO  CALLS  WITH  RFCOMPUTFO  SCALFS. 

JS  T  APT  =  .IP  I  AS  ♦  1 

JENO  =  MINOFI  JH I A9  ♦  J1  ONG «  Jl  AS  T  ♦  JY|  FNGTH  1 
IF (  JSHIFT  .FO.  0  )  GO  TO  ?010 
IF  I  IFxTFNMP)  ,GT.  0  >  GO  TO  POOS 
?00|  IF (  NPAGF  .GT.  1  )  JLTFT  =  0 

IF  (  NPAGF  .FO.  1  )  JLFFT  =  ? 

I F (  NPAGF  .NF.  NPAGFS  )  JR | GHT  s  0 

I F (  NPAGF  .FQ.  NPAGFS  )  JR  I GHT  =  ? 

GO  TO  ?OP0 

Pons  IF  <  ISCALF  .FO.  1  )  GO  TO  pool 

IF<  NPAGF  .FO.  1  .ANO.  IF  XTFNO  ( pi  u.  1  )  JLFFT  *  ? 

I F (  NPAGF  .GT.  1  .OR.  IFXTFNOIP?  .GT.  1  )  JLFFT  =  0 

IF<  NPAGF  .FO.  NPAGFS  .ANO.  IFXTF  '<?)  .FO.  1  )  JR  I  GHT  =  P 

IF  f  NPAGF  .LT.  NPAGFS  .OR.  IFxTFNOI?)  .IT.  3  *  JR  I GHT  =  0 

GO  TO  POPO 

P010  GO  TO  (  P011.P012.P013  )  IFXTFNn(p) 

?0U  JLFFT  =  P 
JR  I GHT  =  1 
GO  TO  POPO 

?0|?  IF (  7SCALFI11  .FQ.  0.00  .OR.  7SCAI F C I )  .PQ.-I.OO  )  GO  TO  2011 
JLFFT  =  0 
JR  I GHT  =  0 
GO  TO  ?opn 
poll  JSHIFT  =  -S 
JLFFT  =  1 
JR  I GHT  =  P 

POPO  IF  (  IF XTFNO 1 1 )  .LT.  ?  >  GO  TO  P100 

IF  (  7SCA|.F<1)  .FO.  0.00  .OP.  7SCAI  F  (  l  )  .F0.-l.00  )  GO  TO  ?100 
IF (  IFX  TFNP  <  ? 1  .LT.  1  .ANO.  NPAGFS  .FO.  1  )  GO  TO  PPOO 

c  this  sfction  gfnfpatfs  thf  hfaofr.  mx  ano  ny  arf  matrix  inofxing 
c  OF  thf  ARRAY  Rt  OCKS  SFmT  SURFP(  OT  im  fxtfndfo  f ALI  s.  npagf  ano 

C  NPAGFS  KCFP  TRACK  OF  THF  PI  OfKs  OF  COLUMNS  DSFO  IN  PRINTING  THF 

C  I  OCAL  ARRAY.  THF  HI  OCKS  OF  C0|  IIMNS  IISFO  IN  PRINTING  THF 


000P3000 
O00P310O 
OOOP3POO 
000P3300 
00023400 
OOOP3SOO 
OOOP3GOO 
00023700 
OOOP3HOO 
OOOP3R10 
OOOP3RPO 
OOOP3B30 
OOOP3B4D 
000P3B50 
O00P3H60 
000P3R00 
00024000 
OOOP4 1 00 
OOOP4POO 
00024300 
00024400 
000P4S00 
OOOP4GOO 
00024700 
OOOP4BOO 
OOOP4QOO 
ooopsoon 
OOOPSIOO 

nnopspno 

000PS300 
0  00PS400 
OOOPSSOO 
OOOPSGOO 
O00PS700 
OOOPSPOO 

ooopsroo 

ooopftoon 

000PG100 

OOOPGPOO 

OOOPGJOO 

ono?G40n 

OOOPGSOO 

nnoPGGOn 

000PG700 

oooPGPon 

OOOPGROO 

OOOPGQIO 

OOOPGRPO 

O00PG930 

000PGR40 


I 


PRINT  9202S.  NX*  NV.  NPAGF  *  NPAGES 
9202S  FORMAT  ( 1HR,  *  (•*  I?  *  *  «*«  I  ?«  *  >  *.  3X  • 

MSURFAfF  PLOT  . PAGE  *♦!?.»  OF  *,I2  > 

GO  TO  2200 

2100  PRINT  92100*  NX*  NY*  NPAGF.  NPAGFS 

921 00  FORMAT ( 1  HR*  //.1H  *  * <*. IP* * . *» I?. * ) • • IX t 

1*SIJRFACF  PLOT  . PAGF  *.12«*  OF  *.IP,  //  > 

PRINT  9?  1 0 1 «  (  ITITLFfl.l.  L  *  1.10  ) 

92101  FORMAT ( 1 H  .  20X.  10<  A«  1,  //  ) 

PRINT  9?10?. 

1  ?  MAX*  MAX  X.  MAX  Y. 

2  7  MTN*  MIN  X.  MIN  Y. 

1  7RANGE.  ITYPF. 

9  7LFVFL  ULFVFLS)  « 

5  7LFVFL  ( 1 )  . 

6  SCALF.  ITYPF. 

7  7  wff, 

fl  7  MIAS. 

9  IXLFNGTH*  JYLFMGTH. 

A  OHSOALF 

92102  F ORM AT ( 1 H  «  20X. 

1 *OA  T A  MAXIMUM  . **Fl  1 .4.1X.*<*.|4,*,*,  14.*)**  /,1H  *  20X, 

2*DATA  MINIMUM  . * *F 1  1  .4  *  1 X  *  * 1 *  *  14 .* .* « 14 .* ) *  *  /,1H  *  20X. 

3*DATA  RANGF  .  *.F11.4*R2,  //*1H  ,  20X. 

4*SCALF  MAXIMUM  .  **F1 1.4*  /,!H  •  20X. 

S*SCALF  MINIMUM  .  *«F11.4*  /»1H  .  20X* 

6*SCALF  RANGE  .  *,Fll.4.R2.  //»1H  «  20X, 

7*Pl.OT  RFFERFNCE  . *.F11,4.  7.1H  .  POX. 

R*PLOT  MT A*.  .  **Fl  1  *9*  7.1H  »  POX. 

9»PL0T  EXTENT  .  (*.I4****, 14**1**  7/*lH  .  POX. 

A*DR  RFFFRFNCE  .  *,FS.2,  //*1H  «  POX. 

R*HA|.F  TONE  DENSITY  WFDGF  AND  SCALF*  ) 

C  THIS  SFCTTON  PRINTS  THF  DFNSITY  WFOGE  OF  GRAY  SHAOFS  ANO  THE 

C  FOllIVAl  FNT  NUMFPICA1  VALUFS*  ROTH  IN  ARSOLUTF  TFRMS  AND 

C  RELATIVE  TO  THE  RFFFRFNC.FS. 

IF  I  LI.F VFI.S  ,GT.  10  )  GO  TO  21  PS 
IF*  ITYPF  .FO.  ORF ILL  )  GO  TO  ? IDS 

no  Pioi  i  =  i.llfvfls 
7LFVFl.lt>  =  7LFVFL1I  1/7MAX 

2101  CONTINUF 

PRINT  9?1 0 1  *  1  7LFVFL1I).  L  =  1  .LLFvFLS  I 
92)01  FORMAT ( 1 H  ,  20X.  J01  2X,FS,3»*MAX*  1  1 

00  2104  l  =1,1  LFVFLS 
7LFVFL1I  )  =  71  FVFL  (L  )*7MAX 
?  1 04  EON  T  Ini  IF 

GO  TO  Plio 


00027000 
000P71 00 
OOOP7POO 
00027100 
00027400 
OOOP7SOO 
00027600 
OO0P770O 
OOOP7ROO 
OOOP7ROO 
00020000 
OOOPrt 1 00 
000P0200 
OOOPHIOO 
00020400 
OOOPOSOO 
OOOPR600 
00020700 
OOOPOflOO 
00020900 
00029000 
00029100 
00029200 
00029300 
00029400 
OOOP9SOO 
00029600 
00029700 
OO029R00 
00029900 
000.10000 
00010100 
00030200 
00010300 
00030400 
00030410 
00030420 
00010430 
00030SOO 
00010600 
00010700 
OOOIOflOO 
000  10900 
000  11000 
00011100 
00031 200 
00031 100 
00011400 

ooonsoo 

00031600 


2G 


?1 0*>  no  210ft  L  =  l.LLEVFtS 

7LEVFLJL)  =  nflSCALF*AL0G10(  7LF VFL <L > /?MAX  1 
21  Oft  CONTINUE 

PRINT  9210ft.  (  7LFVFLCL).  L  =  l.LLFvELS  ) 

9210ft  F ORMAT C 1  H  ,  20X.  101  1T,F7.2«*DR*  >  ) 

00  2107  L  =  1  .LLFVFLS 

ZLFVFLCL)  =  7MAXM  10.00**(  7LFVFL  CL ) /GRSCALF  1  1 
2107  CONTINIIF 

2110  PRINT  92110.  (  ZLFVFL (L) «  L  =  l.LLFvELS  ) 

92110  F ORMAT ( 1 H  .  20X.  IOC  1X.F9.2  >  ) 

DO  21  IS  K  =  1.1 
no  21  IS  L  =  l.LLFVFLS 

IFC  ISOALFCI)  .FO.  0  )  LPRINT  =  CRIAS  ♦  LSTFP*| 

IF  (  ISCALFC3)  .NF.  0  )  LPPINT  =  LRIAS  ♦  LSTFP*<  LIEVFLS»1-L  ) 
JFnGEI  =  21  ♦  C  L- 1  1*10 
JEnr,F2  *  JFOGF1  ♦  9 
no  21  IS  J  =  JFOGF 1 . JFOGF 2 
L INEPLOT (J,K)  =  LINFSVM(LPRTNT.K) 

21  IS  CONTINUE 

no  2iift  j  s  i.2o 

LINFK.M  =  L1NF2CJ)  =  I  INFKJ)  =  RI.ANK 
211ft  CONTINIIF 

JEOGE1  =  JF0GF2  ♦  1 
DO  2117  J  =  JEOGF 1 « 1 3S 
LINEICJ)  =  LINF2IJ)  =  I1NF1CJ1  *  RL ANK 
2117  CONTINIIF 

00  2120  I  =  l.S 
PRINT  92120.  LTNF1 
PR  TNT  02120.  I.INF2 

92120  FORMATClH*.  1  IS (  01  )  » 

PRINT  92171,  LINF1 

92121  FORMAT ( 1 H  ,  11SC  01  )  ) 

2120  CONTINIIF 

GO  TO  ?  J  SO 

2 IPS  00  2 IIS  |.  =  1  .1  LFVFLS 
no  212ft  J  =  1  .  1  IS 

LINFIC1)  =  LINF2CJ)  =  LINF3CJ)  =  RL  ANK 
212ft  CONTINIIF 

IFC  ITyPE  .FO.  OHFILL  1  GO  TO  212R 
7LFVFLCL)  =  71  FVFL  (I  )  /7MAX 
PR  TNT  12127.  7LFVFLCL) 

92127  FORMATC1H*.  21X.  Fft.3.*MAX*  » 

7LFVFLCI.)  =  ZLFVFL  (L  )  *7MAx 
GO  TO  21  30 

2 1  2fl  7LFVFLCI  )  =  0RSCA1  F*AL0G1  0  (  71  FVFL  (L  ) /7MAX  ) 

PRINT  921 ?R.  71 FVFL (L 1 
9212R  FORMATCIH*.  2  I  X  ,F7.2 . *OR»  ) 

7LFVELCI.)  =  7VAXM  10.00**C  7LFVFL  (L )  /ORSCALE  )  ) 

2110  IFC  ISCALFC3)  .FO.  0  )  LPOINT  =  LRIAS  ♦  LSTFP*l 

IFC  ISC  ALF 1 1)  .NF.  0  )  LPOINT  =  LRIAS  ♦  LSTFPM  LLFVFLS* 1-L  1 


00031700 
00031 R00 
00031900 
00032000 
00032100 
00012200 
00032300 
00032400 
OOOlPSOO 
00032ft00 
00012700 
OO032R0O 
00032900 
00033000 
00031100 
00031200 
00033300 
00033400 
OO033S0O 
00033ft00 
00033700 
00033000 
00033900 
00034000 
00014100 
00034200 
00034300 
00034400 
OOO14S0O 
00034ft00 
00014700 
00034000 
00034900 
0OO3S0OO 
00015100 
0003S200 
00015100 
0003S400 
0003SS00 
OOOISftOO 
0001S700 
OOOISROO 
0003S900 
OOOlftOOO 
OOOlftlOO 
0001ft200 
OOOlftlOO 
0001ft400 
OOOlftSOO 
OOOlftftOO 
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<J  u 


no  ?m  k  =  1,3 
no  ?m  j  *  3i,40 

L INFPLOT (  J«K )  =  L INFSYM (LPR1NT»K) 

2131  CONTINUE 

PRINT  42120.  L1NF1 
PRINT  42120*  L1NF2 
PRINT  42121  *  LINF3 
PRINT  4213?.  7LFVEI.IL ) 

4213?  FORMATHH*,  21X.E9.2  » 

PRINT  42120,  L INF  1 
PRINT  42120,  LINF2 
PRINT  42121*  L1NE3 

2I3S  CONTINUE 

this  section  is  nijmrfring  the  columns  un  the  extended  array 

NtlMEROI  OGYI  . 

2150  no  2 1  SI  K  =  1,3 
no  2 ISl  J  *  1 , 1 3S 
L INFPLOT ( J*K  )  s  RLANK 

21  SI  CONTINUE 

JEDGEl  =  JSHIFT  ♦  6 

JEnr,E2  =  JSHIFT  ♦  0  ♦  (  JEND-JSTART  ) 
no  2170  J  =  JSTART »  JFNO 
JPRINT  »  JSHIFT  ♦  7  ♦  J  -  JSTART 
IF  (  J  ,NE ,  MAX  Y  )  GO  TO  21SS 
LINE1 (JPRINT )  =  MEILL 
L INF? ( JPRINT)  =  AE ILL 
LINE3< JPRINT)  =  XEII  L 
GO  TO  2170 

2 1 SS  IF  C  J  ,NE,  WIN  Y  I  GO  TO  2160 
LIUEl  (JPRINT)  *  MEILL 
LINE2(JPRINT)  =  IFILL 
LINE3(  IPRINT)  s  NFILL 
GO  TO  2170 

2160  INDEX  *  J  -  I  J/S  )*S 

IF (  INOFX  .HE.  0  )  GO  TO  2170 
Jinoo  s  (  j/ionn  ) 

Jion  =  (  j-jiooo*iooo  )/ioo 

Jio  =  (  J-Jiono*ionn-Jioo«ioo  1/10 

Ji  =  (  j-jiooo*iooo-jioo*ioo-jio*io  ) 

IF  (  Jiono  .FO.  0  )  GO  TO  2164 

2141  LINE1  (  JPRINT)  =  J100 

2142  L INF? ( JPR TNT )  =  JIO 

2143  L INE3 ( JPRINT )  =  Jl 
GO  TO  2170 

2144  IF (  J100  .EO.  0  )  GO  TO  2165 
GO  TO  2 1 6 1 

216S  L  INF  1  (  JPRINT)  =  DOT 

TF(  JIO  .EO.  0  )  GO  TO  2166 


00036700 

00036000 

00036900 

00037000 

00037100 

00037200 

00037300 

00037400 

00037S00 

00037600 

00037700 

00037000 

00037900 

00030000 

00030010 

00038020 

00030100 

00030200 

00030300 

00030400 

0OO3OSOO 

00030600 

00030700 

00030000 

00030900 

00039000 

00039100 

00039200 

00039300 

00039400 

00039500 

00039600 

00039700 

00039000 

00039900 

00040000 

00040100 

00040200 

00040300 

00040400 

00040500 

00040600 

00040700 

00040000 

00040900 

00041000 

00041100 

00041200 

0004)300 

00041400 
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*  # 


G 0  TO  716? 

7166  L  IMF?  <  IPRIM)  =  HOT 
00  TO  ?163 
7170  CONTINIIF 

PRINT  R7I71.  I  I NF  1 
PRINT  R?l?l.  *  TNF? 

PRINT  R7I71 ,  1  TNF 1 
?  1  RO  no  7  |  9 1  .1  =  JFlTGFl  •  JFOOF? 

L INF1 (  l>  =  MF ILL 
LINE?!  I)  =  WF  ILL 
LINF3II)  =  HLANK 
71R1  CONTINiiF 

IF  (  JIFFT  .OT.  0  )  00  TO  ?1R? 

L  INF  1 (  IFnOFl)  =  HI  ANK 
LINE?!  IFnOFl)  =  HLANK 
?1R?  IF (  JRTOHT  .OT.  ft  )  00  TO  7193 
L INF  1  I  irnOF?)  =  HI  ANK 
L  INF?  I  irnOF?)  =  HI  ANK 
?1  R3  PRINT  R?l?n,  I.  TNF  1 
PRINT  R?l?l,  LINF? 

IF  I  IFXTFNM(l)  .LT.  ?  I  00  TO  ??00 
PRINT  rpirs 
R? I RO  FORMAT I 1  HR  ) 


r  this  is  thf  acihal  priming  sfction  for  the  array  of  fix.yi 

C  VAKIFS.  a  SHARE  OF  gray  IS  ASSIONFO  each  <X*Y)  location  ry 

C  OF  TFRM  T  N I  NO  which  ?|  FVFL  JUST  FXCFFOS  THF  AMPllTlinF  at  that 

C  POINT.  THF  FQIIIVAIFNT  SHAOF  of  GRAY  FOR  THAT  71  fvfl  is  thf 

c  PRINTFO  REPRESENTATION  OF  FIX.Y).  THIS  COMPARISON  IS  OONF  ROW 

c  PY  ROW.  thf  RFMAINOFR  of  this  sfction  is  SPFFOINO  UP  THF  I/O  TO 

C  Thf  l INF  PRINTER, 


??oo  no  ?ioo  i  =  utxifngth 

IX  =  I  ♦  II  AST 

no  ??oi  k  =  l .t 
no  ??m  j  =  i .  ns 

I  INFPl  OT  I  I.K)  3  H|  ANK 
??01  CONTINIIF 

no  ??|S  J  =  JS T  ART . JFNO 
JY  =  J  -  J|  AST 

INOFX  =  l  *  (  JY-I  ) •  l  XI FNGTH 
HO  ??|4  L  =  I .LLFVFI  S 
IF I  L  .FO.  LLFVFI S  )  00  to  ??ns 
7TFST  =  70ATA ( INOF  X  > 

IF  I  ISC  ALE  (4)  .FO.  0  )  GO  TO  ??0? 
7TFST  =  ITFST 

??0?  IF!  7TFGT  .OT.  71  FVFL  II  I  )  00  TO  ??14 
??0S  JPRINT  =  ISHIFT  ♦  7  ♦  J  -  JST ART 


00041 SOO 

00041600 

00041700 

0004 1 POO 

00041900 

00047000 

0004?1 00 

0004??00 

00047300 

00047400 

0004?S00 

0004P600 

00047700 

00047FIOO 

00047ROO 

00043000 

00043100 

00043700 

00043300 

00043400 

00043600 

00043600 

00043700 

00043R00 

00043900 

00043R1 0 

00043970 

00043930 

00043940 

00043950 

00043960 

00043970 

00044000 

00044100 

00044700 

00044300 

00044400 

00044600 

00044600 

00044700 

OOO44AOO 

00044900 

00045000 

00046100 

0004S700 

0004S300 

00046400 

00046600 

00045600 

00045700 


29 


IF (  fSr«LF(l)  .FO.  0  »  LPRINT  =  LRIAS  ♦  LSTFP*l 

IF(  IsrALFIl)  .NF.  0  1  LPRINT  =  LRIAS  ♦  LSTFP*!  LLFVELS*1-L  ) 

no  ?21  0  K  =  1.1 

L  INFPLOT (  (PRINT. K>  =  L INFSYM 1LPR INT ,K 1 
2210  CONTINUF 

00  TO  2215 

2214  CONTINUF 

2215  CONTINUE 

222 0  IF  (  JLFFT  .FO.  0  1  00  TO  2221 
LINF1  (JFOGFl )  =  MF1I.L 

L  INF?  < JFOGF  1 )  =  WF  ILL 

2221  IF (  JRIOHT  .FO.  0  )  GO  TO  2230 
LINF1  I JFOGF?)  =  MF  ILL 
LINE2UFDGF2)  =  WF  ILL 

2230  IF  I  IX  .NF.  MAX  X  )  GO  TO  2240 
IF  (  JLFFT  .LT.  2  )  GO  TO  2231 
L  INF  I  CJSHTFTM)  *  MFILt 
LINF1 ( JSHIFT ♦  ?)  *  AFILL 

I  INF  1 (  ISHIFT*3>  s  XFILL 

??31  IF (  JPTGHT  .LT.  ?  >  GO  TO  2232 
L INF  1 ( JF0GF2*3)  *  MF ILL 
LINF1 <JEDGF2«4)  =  AFILL 
LINF1  (  IF0GF?«5>  =  XFILL 
GO  TO  2260 

??32  IF (  L INFl ( JSH IFT ♦ 1 )  .FO.  MFILl  )  GO  TO  2260 

2240  IF  (  IX  .NF.  MIN  X  )  GO  TO  2250 
IF  (  JLFFT  .LT.  2  1  GO  TO  2241 

L  INF  1 (JSH IFT ♦ 1 )  =  MFILL 
L INF  I I JSH IFT ♦  ?)  s  IF  ILL 
L INFl (JSH IFT ♦ 3 )  =  NF  ILL 

2241  IF  C  JRTGHT  .LT.  2  )  GO  TO  2242 
LINF1 (JF0GF2*3)  *  MFILL 

LINF1 (JFDGF2*4)  =  IFILL 
LINEl (JFDGF2*5)  *  NF ILL 
GO  TO  2200 

2242  IF (  L INF  1 ( JSH I FT ♦ 1 )  .FO.  MFILL  )  GO  TO  2260 

2250  IF  C  JLFFT  .LT.  2  .ANO.  JRIGHT  .LT.  2  1  GO  TO  2260 
INDEX  =  IX  -  (  IX/5  )«S 

IF (  INOFX  .NF.  0  )  GO  TO  2260 
11000  *  (  T X/ 1000  ) 

1100  *  <  IX-I1000*1000  1/100 

110  =  (  TX-I1000*1000-I100*100  1/10 

II  =  (  IX-I1000*1000-I100*100-I10*|0  1 

IF (  JLFFT  .EO.  2  )  ASSIGN  2251  TO  JNEXT 
IF  C  JLFFT  .LT.  2  1  ASSIGN  2250  TO  JNEXT 
GO  TO  (NEXT,  <  2251 .2259, 2260  1 

2251  INDFX  3  JSHIFT  ♦  1 
ASSIGN  2259  TO  JNFXT 

225?  IF (  11000  .FQ.  0  )  GO  TO  2256 

l infi (Index i  *  nooo 


00045800 

00045900 

00046000 

00046100 

00046200 

00046300 

00046400 

00046500 

00046600 

00046700 

00046800 

00046900 

00047000 

00047100 

00047200 

00047300 

00047400 

00047500 

00047600 

00047700 

00047B00 

00047900 

00048000 

00048100 

00048200 

00048300 

00048400 

00048500 

00048600 

00048700 

00048800 

00048900 

00049000 

00049100 

00049200 

00049300 

00049400 

00049500 

00049600 

00049700 

00049800 

00049900 

00050000 

00050100 

00050200 

00050300 

00050400 

00050500 

00050600 

00050700 
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??5i 

l.  INF1  IlNnFXM  )  =  MOO 

OOOSOROO 

PPS4 

LINFi  (!Mnrx«?)  =  no 

00OSO900 

PP55 

L  INF  I  I  TNnFx ♦ 1)  =  11 

ooosiooo 

GO  TO  IMF  XT  *  1  ??S1 ,??S9,??40  ) 

OOOS1 100 

PPS4 

l  I NF  1  (  T Mt )F  X  >  =  OOT 

ooosi?oo 

IF  (  lino  .FO.  0  )  GO  TO  ??57 

00051300 

GO  TO  ??S1 

00051400 

??5  7 

L INF  1 ( INOFX ♦ 1 1  =  OOT 

0005)500 

I F  <  11  n  .FO.  0  1  GO  TO  PPSR 

000S1400 

GO  TO  ?P54 

00051700 

??SR 

L  INF  1 1 I MOF X  ♦? )  =  OOT 

00051 ROO 

GO  TO  ??S5 

00051900 

??5Q 

I F (  JR1GHT  .LT.  ?  )  GO  TO  ??40 

0005?000 

ASSIGN  ??40  TO  INF  XT 

0005P1 00 

INOFX  =  JFOGFP  ♦  ? 

0005??00 

GO  TO  ??S? 

0005P300 

??40 

PO  ??70  K  =  1,1 

0005P400 

no  ??Tn  J  =  1.17 

0005P500 

I F (  J  .FO.  1  1  GO  TO  ??45 

0005P400 

INOFX  =  I  J-?  ) *H  ♦  7 

0005?700 

LINFP1.0T(IM0EX*1,K)  =  ISHTFTl*LTNFPI.OT  (INOFX*l»K) 

0005?fl00 

LINFPLOT I INOFX*?. K>  =  1  SH  IF  T  ?*L  INFPI.  0  T  ( INOFX *?,K) 

0005P900 

l  INFPIOT  (  TNOFX*1,K)  =  1  SHIT  T  1*1  INFPI  OT  t  INOFX ♦  1  .K ) 

00053000 

L INFP!  OT< TNOFX*4,K)  =  I5HIF  T 4 *L INFPLOT (  INOFX  *4, K ) 

00053100 

L INFPIOT ( 1MOFX.S.K)  r  I SH IF T S«L INFPI OT  1  INOFX *5 *K ) 

noosipoo 

L  INF  PLOT  (  !NnFX*4»K)  r  ISHIFT4*L  INFPI  OT  C  IN0FX*4.K) 

00051300 

1  INFPIOT (TM0FX*7,KI  =  I SMIF T 7*L INFPI  OT ( I NOFX ♦ 7 .K ) 

00053400 

JLlNF  (  (,K  1  =  1  INFPlOT(INOFX*1  ,K>  .OR.  1.  INFPLOT  ( INOFX*?, K> 

00053500 

1 

.OR.  L INFPI  OT ( INOFX ♦  1  «K )  .OR.  1  INFPLOT ( INOFX *4, K  ) 

00053400 

? 

'  .OR.  LINFPLOT  (I  NOFX  *5  »  K  >  .OR.  1.  INFPL  OT  ( INOFX*4»K) 

00051700 

1 

1  .OR.  L INFPI  OT C1N0FX*7.K>  .OR.  1  INFPLOT ( INOFX  *R»K ) 

OOOS3ROO 

GO  TO  ??70 

00053900 

PP45 

IF  (  K  .FO.  3  )  IPAfiF  =  ISHIFT I *RLANK 

00054000 

IF  (  K  ,MF.  1  )  IPAGF  =  ISHIFT1  *PLIIS 

00054100 

LINFPLOTII  .*>  =  ISHIFTPHINFPI  OTI1.KJ 

OOOS4POO 

1  INFPLOT!?, K)  =  ISHTFT1*!  INFPLOT!?, K) 

00054300 

L  INFPLOT  !  1,10  =  ISH1F  T 4*1.  INFPI  OT  ( 1,K  > 

00054400 

L  INFPLOT  14. K)  =  ISHIF  TS*L  INFPI  OT  14. K  ) 

00054S00 

L  INFPI  OT  IS.K)  =  ISHIFTG«|.  INFPLOT  IS. K) 

00054400 

L INFPLOT 14. K)  r  ISHIFT7»LINFP|  OT 14, K> 

00054700 

JLINFIJ.K)  =  IPAGF  .OR.  LINFPLOT  II ,K) 

OOOS4ROO 

1 

.OR.  LINFPI  OT !?,K)  .OR.  L INFPLOT  1 3, K » 

00054900 

P 

.OR.  L INFPI OT 14, K)  .OR.  L INFPLOT  IS, K ) 

00055000 

3 

i  .OR.  1  INFPI OT 14, K)  .OR.  L INFPLOT ( 7,K) 

00055)00 

??7n 

CONT  INOF 

0005SPOO 

PRINT  9??70,  J|  INF 

00055100 

9PP70 

FORMATI  171  AR  )  ) 

000SS400 

P100 

CONTINMF 

00055500 

00055400 

C 

THIS  SrrTTON  IS  NIIMHFRING  THF  COlllMMS  I  IN  FXTFNOFP  ARRAY 

00055410 
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c  NUMFROLOGY  J  . 

DO  2101  K  =  1.1 
DO  2301  J  =  1.135 
L1NFPIOT (J.K)  =  BLANK 
2301  CONTINIIF 

IF (  IFXTFNO ( 1 1  .FO.  0  .OR.  IFxTFNO < 1 1  »FQ.  3  1  JTAIL  «  1 
IF (  IFXTEND < 1 1  .FQ.  1  .OP.  IFXTFNO(l)  .FQ.  2  )  JTAIL  «  0 
IF  (  JTAIL  .EQ.  0  )  f,0  TO  2400 
JEOGF1  =  'SHIFT  ♦  6 

JEDGF2  =  JSHIFT  ♦  8  ♦  <  JFNO- JSTART  ) 

00  2202  J  =  JFOGE I  ♦  JEDGE2 
L INF  1 (J)  =  MF ILL 
LINF2U)  =  WFILL 
230?  CONTINHF 

IF  (  JLFFT  .GT.  0  )  GO  TO  2303 
L INE 1 ( JFOGF 1 )  =  BLANK 
L  INF2 ( JFOGF 1 )  =  BLANK 

2303  IF  (  JRIGHT  .GT.  0  )  GO  TO  2304 
L  INF  1 (JFOGF?)  s  BLANK 

l INF? (  IFOGF2)  r  BLANK 

2304  PRINT  92120,  LTNF1 
PRINT  02121,  LINE? 

L INE 1 ( JFOGF 1 )  =  LINF1 (JEOGE2)  =  BLANK 
L  INE2 (JFOGF 1 >  =  L !NF?( JFOGF2)  =  BLANK 
00  2320  J  =  JSTART.JFNO 
JPRINT  =  JSHIFT  ♦  7  ♦  J  -  JST ART 
LINF1  (  (PRINT)  s  LINF?(JPRINT)  =  BLANK 
IF (  J  .NF.  MAX  Y  )  GO  TO  2305 
LINE) (JPRINT)  =  MF ILL 
LINF?(JPR!NT>  =  AFIL) 

LINF3<  (PRINT)  s  XFII  1 
GO  TO  2320 

2305  IF (  J  .NF.  MIN  Y  >  GO  TO  2310 
LINF1 (JPRINT)  =  MFILL 
LINF2( JPRINT)  =  IFILl 

L INF3 ( (PRINT)  =  NF ILL 
GO  TO  2320 

2310  INOFX  =  j  -  (  J/S  )*5 

IF  (  INOFX  .NF.  0  )  GO  TO  2320 
J1000  5  J/1000 

J100  =  (  J-J1000*1000  )/100 

J10  =  (  J-J1000*1000-J100«100  1/10 

J1  =  (  J-J)000»1000-J100*100-J10*10  ) 

IF (  J1000  ,E0.  0  )  GO  TO  2314 

2311  LINF) ( IPRTNT)  =  J100 

2312  I  INF? (JPRINT )  =  J10 

2313  LINF3( (PRINT)  =  J1 
GO  TO  2320 

2314  IF (  J100  .FO.  0  )  GO  TO  2315 


00055620 

00055700 

00055800 

00055900 

00056000 

00056100 

00056200 

00056300 

00056400 

00056500 

00056600 

00056700 

00056800 

00056900 

00057000 

00057100 

00057200 

00057300 

00057400 

00057500 

00057600 

00057700 

00057800 

00057900 

00058000 

00058100 

00058200 

00058300 

00058400 

00058500 

00058600 

00058700 

00058800 

00058900 

00059000 

00059100 

00059200 

00059300 

00059400 

00059500 

00059600 

00059700 

00059800 

00059900 

00060000 

00060100 

00060200 

00060300 

00060400 

00060500 
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o  r> 


V  / 


00  TO  ?111 

?31S  L INF  1 (  JPPTNI)  =  OOT 

IF (  Jin  .FO.  0  )  00  TO  ?316 
00  TO  ?11? 

?116  l  INF?( JPR1NT1  =  OOT 
GO  TO  ?111 
?3?0  CONTINUF 

PRINT  R?l?l .  t INF1 
PRINT  R?l?l ,  t INF? 

PRINT  R?l?l.  L  INF1 

?400  JR  IAS  =  JFND 


?son  continof 

THIS  FNr>S  THF  PRINTING  SECTION. 

THF  FXTFNOFO  ARRAY  PARAMFTFRS  ARF  IIPDATFD  IN  THIS  SECTION. 


IF  (  !FXTENnM>*!EXTFNr><?)  .FO.  0  )  00  TO  30?0 
IF  I  lEXTFNn<l)«IEXTFNnl?»  .FO.  <J  )  00  TO  1015 
00  TO  (  lool tinns.ioio  >  toxy  ♦  l 

inm  IlAST  =  IX|  FNOTH 
JL AST  =  JYl FNOTH 
NX  =  0 
NY  =  0 
00  TO  SOOO 

ions  IF  (  IFXTFNnm  .FO.  1  )  00  TO  1000 
ILAST  =  ILAST  ♦  I  XL  FNOTH 
NX  =  NX  ♦  1 

GO  TO  5000 
3006  ILAST  =  0 

JL AST  =  JL  AST  ♦  JYl  FNOTH 
NX  =  1 
NY  =  NY  ♦  1 
00  TO  1010 

1010  IF (  IPxTENn (?)  .FO.  1  1  00  TO  3011 
JL  AST  =  JL AST  ♦  JYl  FNOTH 

NY  a  NY  ♦  1 
00  TO  SOOO 

1011  JL AST  =  0 

ILAST  =  ILAST  ♦  I  XL  FNOTH 
NY  =  I 
NX  =  NX  ♦  1 
GO  TO  1015 

1015  ILAST  =  JL AST  =  0 
NX  =  NY  =  1 
IOXY  =  0 

3015  PRINT  RIO  1  5 


00060500 
00060700 
00060000 
00060900 
00061000 
00061100 
00061?00 
00061300 
00061400 
00061500 
00061600 
00061700 
00061  POO 
00061 ROO 
0006?000 
0006P100 
0006P1 10 
0006?1 20 
0006??00 
0006P300 
0006?400 
0006?500 
0006?600 
0006?700 
0006P800 
0006?R00 
00063000 
00061100 
00063P00 
00063300 
00061400 
00063500 
00063600 
00063700 
00061800 
00063ROO 
00064000 
00064100 
00064?00 
00064300 
00064400 
00064500 
00064600 
00064700 
00064800 
00064R00 
00065000 
00065100 
00065?00 
00065100 
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93016  FORMAT  ( |HA»  //  1 
50  TO  5000 

30?0  IF  (  IFXTFNOm  .NF.  0  1  50  TO  3040 
IF  <  IFXTFNO (?)  .NF.  0  >  GO  TO  3030 
GO  TO  3015 

3030  IF  t  IFXTEnM?)  .EO.  3  )  GO  TO  3015 
JLAST  *  Jl AST  ♦  JYLFNGTH 
NY  a  NY  ♦  1 
GO  TO  5000 

3040  IF  1  IFXTENOl 1 )  .EO.  3  )  GO  TO  3015 
IL»5T  =  HAST  ♦  I XLENGTH 
NX  a  NX  ♦  1 


5000  RETURN 


FNO 


00065400 

00065500 

00065600 

00065700 

00065500 

00065ROO 

00066000 

00066100 

00066?00 

00066300 

00066400 

00066500 

00066600 

00066700 

00066500 

00066900 

00067000 

00067100 

OO067POO 

00067300 

00067400 
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COMPARISON 


No  equivalent  plotting  routine  was  found  in  the  library. 


7.0  TEST  METHOD  AND  RESULTS 


The  following  simple  example,  PR0GRAM  TEST,  illustrates 
the  use  of  SURFPL0T . 


PROGRAM  TFST 

C  THTS  IS  AN  FXAMPLF  PROGRAM  USING  SURF  PLOT.  THF  FUNCTION  BF1NG  PLOTTER  IS  THE 
C  SUM  op  TWO  SINUSOIDS.  ONE  IN  X  ANO  ONE  IN  Y.  ON  A  HIAS. 

OIMFNSTON  A (77, 1 00) 

OIMFNSTON  70ATA  (7?. 100)  .  ITITLEIIO),  ISCALE  (4)  »  7SCALEO) 

DIMENSION  TEXTFND  <7) «  7LFVEL (10)*  7PARAMTR (6) 

1XLFMGTH  =  7? 

JYLFNGTH  =  100 

READ  90000,  (  ITITLFm*  1  =  1.10  ) 

90000  FORMAT (  10<  AA  )  ) 

PI  =  1. 141S9PGS1S 
RADIAN  X  =  PT/17.0 
PARIAN  Y  *  PI/70.0 
DO  1000  J  =  1,100 
ARG  Y  =  J*RAD1 AN  Y 
PART  Y  =  S1NF (  ARG  Y  ) 

DO  1000  I  a  1,77 
ARG  X  a  T»RA01AN  X 
PART  X  =  SINF (  ARG  X  ) 

AII.J)  =  7.00  ♦  PART  X  ♦  PART  Y 
1000  CONTINUF 

CALL  SURF  PLOT  (  A,  I X«  FNGTH,  JYLFNGTH,  ITITLF •  ISCALE.  7SCALE ♦ 

1  7LFVFL.  7PAPAMTR.  IFXTFND,  ISFARCH  ) 


FND 


THF  FUNCTION  TS  F(X,Y)  =  7.00  ♦  S INF <PI *  I / 1 7 . 00 )  ♦  S INF (P I  * J/70 . 00 > 
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<  1«  li  Surface  plot 


PAGE  1  OF  1 


THE  FUNCTION  IS  F|X,t|  ■  2,00  ♦  SINE  |  P  !•  1/ 12.00  >  •  SINF  |PI*  J/20 ,00) 


DATA  MAXIMUM 
OATA  MINIMUM  .... 
OATA  *ANGC  . 


SOLE  maximum  .... 
SCALE  minimum  •... 
SCALE  range  ••«... 


PLOT  REFERENCE 
PLOT  HAS  ..... 
PLOT  EXTENT  ... 


OB  REFERENCE 


4.0000*000  I 
0.0000*000  ( 
♦. 0000*000 

4.0000*000 

4.0000-001 

3.6000*000 

4.0000*000 
0.0000*000 
l  ?2*  100) 


6  • 
l  8  • 


10) 

30) 


0.00 


half  tone  density  -EOGE  AND  SCALE 

O.lOOMAa  o,t  UMAX  O.JoOMAx  0.A00MAX  O.^OOHXX  0.600MAX  O.TOOMAX  0.800***  O.OOOMAx  l.OOOMAX 

4.00-001  0.0-  -001  1.20*000  1.60*000  2*00*000  2.*0*000  2.80*000  3.20*000  3.60*000  4.00*000 

... - ••••••••••xxxxxxxxxxooOOOOOOOOBatMMMOtMIBtaiBBMVVlBVVMlMHHtM 

•*-———♦•♦•♦•♦•  ♦♦••••••••••XAXXXXXXXxOOOOOOOOOO«MMMMiamtMM6«MV«9fVt«*IM>HMt 

. . .  ♦  •  .*•••••••  •••xixxxx  xx  a  AOOOOOOOOuUMBItMMaitf|aaaa>Mtl|fffMBIMMtBBBB 

—  •  —  •*.*••  ••••••••••••••xx  xxxxxxxxOOOOOOOOOUBtSaBMaiMItMtlMMHMVVMMBMtMMt 

.••...•—••«•••••  ..••••••••••xxxxxx  xxx  XOOOOOOOOOOIIIIinBINNHHIIMNIHIMlINNNNI 

•  M  ,  •  .  M  ,  ......1 

.  A  1  2  2  I  3  4  4  9  S  6  6  T  7  0  A  0  R  o 

5X5  05N5050S05  OSOSr  50 


••15 

MIN 

..20 


.25 


■•••■•t«V0«V|IIIMMoOXAX4444««*»«««xxxoOR|0Mt0fMM0MIIMOOXXA**44»«4444MXXOO«tMIMV«MMM0M«OM 

M10«VVV90MMH0aa«0OOOXXXX4«44*xiXXOOOM0«M9VI«l00VVVaMIIOOOAXXX»«*»*xxxXOOOMM09Hm«MfVMM0N 

NaiaaaaaMHavaavaaaaaouooxxKxxxxxxxoooaiaaaaaaaaMaaaaaaaaMaooooxxxxxxxxxxoooaaaaaavMHMaaaaaaBaaa 

!!aaH>*tii»uaBaaaaaaaaaaoouooxuooooaaaaaaa!aaaaiaataaaaaa«!aaaaaoooooxoooooaaaaaaaHiiaH»SMMaa«aa!! 

NaaaaaaaauataaaaaaBaaaaoooooxxxoooooaaaaaaaamaaaHaaaafvaaaaaoooooxxxoooooaaaaaaaaiMaaHNMaaaaaa 

RaaaaaB»MH»Ma«aaaaaaaoooxKxxAxxoooaaaaaaB«aiaBaHMiai«aaaaaaoooxxxxxxxoooaaaaaa«aiHmmHaaaaaH 

MaaaaavailHaaaavaaaaaooooxxxxxxxxxxooaaaaaaBaaflHUaBavaaaaaooooxxxxxxxxixoooaaaaaaaaviBiMaiaaaaaaH 

Maaaif«fV««vvviaaaaaoooxxxx444*4xxxxoooaoaaa9aavvavvafaaaaaaoooxxxx«**«4xxxioooaaaaa99avaaaB«aaaaaaaH 

Haaaaaaafi99«aaaaaaaooxxx4««»*««444«xxxooaBaaaaaava«vaaaaaaoooxxx444**44«»»«Kxxooaaaaa8a9aa«vaaaaaaaoM 

Mooaiaaaaaaaaaaaaaooo8x«4*****.*.***«4xxxooaaaaaaaaaaaaaaaoooxx«************44xxxooaa»aaaaaaaaaaaaoooH 

"xoooaaaaaaaaaaaoooxx**** . — •••••••ixouoaaaaaaaaaBaoooxx*******- . .♦•••xx0ooa»i»aaaaaiat)oo*ia 

■xxxooooaaaaoooooxxx*****— ——— ---♦♦♦••xxxooooaaia»ooooxxx##*  ♦••xx  xoooosaaaaoooox  xx#" 


•  *5 

MAX 


R4XXXXOOOOOOOOOOX 

■•••xxx0000000*xxp****»“— 

N« • •XXXXXOOOX XX AK**« 

■••••xxxxxxxxxxx***..*— 

■  •••XXXXXOOOXSXXX***.*  — 

R***xxxoOOOOOO*xx*##***-*- 

"•ixxxOOOOOOOOOxxxx*****--— 

■xxxooooaaaaaooooxxx*****— 

axoooaaaaaaaaaaaooos 


••xAXXOOOOOOOOOO* 

..... ♦  •••XXXOOOOOOOX  A  X»****  —  - - 
•.••.•••XXXXXOOOAXXXX••••••-- 

•-•••••••xxxxxxxxxxx******--- 

---..•••XXXXXQOCXK XI ■••*♦•-- - 
-..♦.♦•••XXXOOOOOOOXXX******—- 

--••♦*• ••xxxxooooooooooxxx*****-- 
- - ••4xxxooooaaaaaooooxxx**.  ♦  ♦- 

xoooaaaaaaaaaaaoooxx**** 


ooooooooooxxx«»a 
--••  •  ♦  •  •  •  x  xxooooooox  X  x  •  ••  .  ii 
---♦♦♦••XXXXXO  00  XXXXX****M 
---•♦♦•••XXXXXXXXXXX4*4.t|( 
— -• • • ••xxxxxooox  X  XXX****M 
— ♦♦••••XXXOOOOOOOXXX»4«*K 
.  .... ••♦••XXXXOOOOOOOOOOXXX**" 

------- •••••*xxooooaaaaaooooxxx*n 

♦— . ♦••••xxoooMaaaaaaaaaoooxxM 

HOoaBaaaaaaaaaaaaaoooRx««4******«***«4xxxooaaaaaaaaaaaaataoooxx**4»*** . •••xxxooaaaaaaaaaaaaaaaooo" 

. . . 

"aaaaaB«aaa«viaaaaaaaoooxxxx444P*xxxxoooaaaaaaa«a«avtaa«aaaaioooxxxx»4444xxxxoooaaaaaa«vv«iavvaaaaaaa" 

HaaaavnHmavvvaaaaaoQOoxxxxxxxxxxuouaaaaaBaftBMaiaaaaaaaaaooooxxxxxxxxxxoooaaaaaavviHlBBaaaaaaaBM 

"aai«aMUimi«fvaaaaaaoooxKxxxxxoooaaaaaaavvM»aaataBVfiaaaaaaoooxxxxxxxoooaaaaaaaaaiiHauB»anaaa" 

"aaaatHMMiMtavaaaaaaoooooxxxoooooaaaaaaaaaBaaaaaaBBaaaaaaaaaoooooxxxoooooaaaaaavaBMBaBBMMaaaaaM 


■aaaBaaaaMaiaBaaafiaaaaoooooxxxoooooaaiaaaataaaaBaaaBBaaaaaaaiaoooooxxxoooooaaaaaaaBaBiHHMBaaaaaa" 
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Haaaaaaaavvfiaavaaaaaoooxxxx*»*«»xxxxoooaaaaa«aaa«avvafvaaaa»oooxxxK*«*««xKKxoooaaaaaaaaiaaa«aaaaaaaa" 
"aiaaaaaav9fiMtMaaooxxx4**4*4*4«««xxxooaaaaaaa«««vvaaaaaa0ooxxx4*»»**««*4«xxxooaaaaaaamaaaaaaaBao" 

. . ..••••xiooooaaaaoaaiaaaoaaoooii********* . xxxooaaaaaaaaaaaaaaaooo" 

«xoooaaaaaaaaaaaoooix#******— — ♦••••••xxouoaaaaaaaaaaaoooxx«»»*«  ♦♦—...♦•  •••••xxoooaaaaaaaaaaaoooxxM 

"xxxooooaaaaaaoooxxx***  . . — — — ♦••••xxxcoooaaaaaooooxxx4ii 

"•XXXOOOOOOOOOOOXXX*****— •  —  —  ••  ••XAXOOOOOOOOOOOXXX***  — —  —  —  ••♦••XXXXOOOOOOOOOOXXX**" 

"•••IXXOOOOOOPXXX«»«***-—  — .  ... •••XXXOOOOOOOX  XX«*»***—  —  •• ••••XXXOOOOOOOX XX*««*M 

"•••xxxxxoooxxxxx***  •  •---  ... ..•••xxxxxoooxxxxx*****— •  ... •••••xxxixoooxxxxx****M 
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"•••xxxxxoooxxxxx***  •  ♦•••  - . •XXXXXOOOX  . . 

^•••XXXOOOOOOOXXXP**#**—  — .  ♦•♦•••xxxoooooooxxx******  — 

•  •45  «*XXXX000000000XX1XP»***— —  — ' —  •••••xxxxoooooooooxxxx***..— 

. . — — -♦•♦••xxxooooaaaaaooooxxxt*.*.- 

"xoooaaaaaaaaaaaoooxx*******— — •••••••xxoooaaaaaaaaaaaoooxx****** 
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. . . . •••••xxxooaaaaaaaaaaaaaaaooo" 

"aaaa8aaa9aaaaaaaaaaooxxx««**44**«««xxxooaaaaaaaaaaaaaaa8aaaooxxx»*4»4»44444xxxooaaaaaaaaaNaaaaaaaao« 

•  •so  »aaaaaBfV99a«if9aaaaaoooxKxx*«*««xxxxoooaaaaat«9a9aaaaaaaaaaaooo*xxx»»*44xxxxoooaaaaaaaaiaaaavvaaaaja" 

aaaa  vavaaaaaaaaa  aaaiaaoooo  x  x  x  x  x  x  x  x  x  xoooaaaaaaaaaaaaaaaaaaaaaaaooo  ox  x  x  x  x  x  x  x  x  x  oooaaaaaaaaaaaaaaaaaaaaaa" 
"aaaavMMHaMfaBaaaaaaoooxxxxxxxoooaaaaaavvaMHaaBaBaBaaaaBiaoooxxxxxxxoooaaaaaavaaaaaaaatMaaaaaaM 

"aaaoaaaaaaaBiaaiiaiaiaatoooooxooooooaoaaafvaaaaaaaaaaaaaBiaaataaoooooxoooooaaaaaaaaaaaaaaaaaaaaaaaaa" 

•  •55  "avaattMaHMHaavaaaaaoooooxKxoooooaaaaaftaMaaatMaaaivvaaaaaoooooxxxoooooaaaaaaaaaaaaaNNMaafaaM 

HaaaaiHaiMaaMavaaaaaaoooxxxxKxxoooaaaaaaaaaMBBaaaaaviaaaaaaaoooxxxxxxxoooaaaaaaaaaNHamBfBaaaa" 

"aMaaaaatNaaaaaaaaaaooooxxxxxxxxxxoooaaaaaataaaaaaaaaaaaaaaaooooxxxxxxxxiioooaaaaaaaaaHatBaaaaaaaa" 

"aaaaaa«vfaav««aaaaaaoooxxxx»444«xxxxoooaaaaaaaaaaaaiaaaaa8oaoooxxxx»»**«xxxxoooaaoaaaaaaaaaaaaaaaaaa" 

"aaa8aaaaaoaaaaaaiaaooxxA«44»4*4«4«4xxxooaaaaaaaaaaaaaaaaaaaooxxx«**««»*44*«xxxooaaaaaaaaaaaaaaaaaaao" 

•  •60  . . . . •••xxxooaaaaaaaaaaaaaaaooo" 

"xoooBuBV;MBBBaoooix««*****— ^ - . •••xxoooaaaaaaaaaaaoooxx*##.*.*— — . •xxoooaaaaaaaaaaaoooxxM 

"xxxoooocsaaaooooxxx««***— — — — . . . — — — . ••xxxooooaaaaaooooxxx«M 

"•XX XOOO.OOOO0OOXXXRM*  ••----  ----•♦•••XXXXOOOOOOOOOOXXX*****  — —  .... •••••XXXXOOOOOOOOOOXXK**M 

M»»*XXXOOOOOOOXXX««»***—  ...  ••••••XXXOOOOOOOXXX******—  . .  XXXOOOOOOOX  XX****M 

•  •65  M*4*XXXXX0*)0XXXXX»#* * •-•-  . . .  . •••XXXXXOOOXXXXX»«.*M 

"♦•••XXX. ’X  XXXXXX ••••♦• —  —  —  •— •••••••XXX XXXXX X XX******  — •  ---•♦••••XXXXXXXXXXX**4**M 

"♦••XXXXXOOOXXXXX*****  — -  ... •••••XXXXXOOOXIXXR*****—  —♦•*••  X  X  X  X  XOOOXX  XX  X****M 

■•••XXXOOO'.OOOXXX***#**—  ...  ••♦•••XXXOOOOOOOX  xx»*****  — -  —  •♦♦••*XXXOOOOOOOXXX»«**M 

M»XXXXOOCO;OOOOXXXX*****  —  —  —  —  •••••XXXXOOOOOOOOOOXXX  •••♦•  — --  —  —  •••**XXXXOOOOOOOOOOXXX**M 

•  •to  "xxxooooaaaaaaoooxxx*****— — — . . . ------ . . 

. . .  •••••iioooaaaaaiiiaiaonox*****  *♦.----- •••••••xxoooaaaaaaaaaaaoooxxM 

"ooaaaaaaBiaaaaaaaooxxx«*««******.«***xxxooaaaaaaaaaaaaaaaooxKX4««*«**«**«««4*xxxooaaaaaaaaaaaaaaaooxM 
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8 . 0  REMARKS 


The  ''writing  speed"  of  the  line  printer  is  not  high. 

The  technique  used  here  depends  strongly  on  the 
natural  spatial  integration  done  by  the  eye.  If  the 
structures  change  quickly,  e.g.,  in  the  span  of  one 
element  to  its  neighbor,  then  the  detail  will  not  be 
perceived. 

The  gray  scale  used  in  SURFPL0T  is  one  selected  after 
some  trial  and  error.  It  is  by  no  means  unique.  There 
is  literature  on  the  subject,  although  the  only 
reference  consulted  was  Knuth  [1].  The  user  who  wants 
to  develop  a  different  gray  scale  can  either  (1)  use 
more  than  ten  shades  in  which  case  this  subroutine 
requires  a  new  gray  scale  or  (2)  change  the  internal 
gray  scale  contained  in  the  data  statements  for  LINESYM. 
One  technique  of  interest  is  to  "go  around  the  gray 
scale"  several  times  to  obtain  more  shades.  That  is, 
using  the  same  symbols  in  cycles.  The  author  would  be 
interested  to  see  other  gray  scales  generated  using 
the  line  printer. 
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